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ABSTRACT 


Evidence is presented for the X-ray variability of Cyg X-l on time 
scales down to a millisecond. Several "bursts" of millisecond 
duration are observed. The duty cycle for bursting is estimated to be 
^2 x 10 averaged over the entire 49. second exposure, although the 
maximum burst activity is associated with a region of enhanced emission 
lasting about l/3 second. Such bursts may be associated with turbulence 
in disk accretion at the innermost orbits for a black hole. 
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MILLISECOND TEMPORAL STRUCTURE IN CYG X-l 
I. INTRODUCTION 

Millisecond and/or submillisecond pulses from an X-ray source have 
been advanced as possible observational evidence of a black hole (Shakura 
and Sunyaev, 1973, Leach and Ruff ini, 1973, Pringle and Rees, 1972). 

Cyg X-l is already known to display temporal structure on time scales down 
to several tens of milliseconds (Oda, et al., 1971, Rappaport, et al., 

1971, Holt, et al., 1971, Schreier, et al., 1971). In this Letter we 
report the observation of millisecond bursts from Cyg X-l. 

II. EXPERIMENT 

On October 4, 1973 at 340 UT an X-ray astronomy payload was launched 
(flight number 13.010) from White Sands Missile Range in order to study 
the compact binary sources Her X-l, Cyg X-l and Cyg X-3. The payload 
consisted of two multilayer multianode gas proportional counters with 
total effective viewing area of 1360 cm and with 1-mil aluminized Kapton 
windows. Both detectors were filled to 1.1 atmospheres, one with 90% 
argon - 10$ methane and the other with 90$ xenon - 10$ methane. The 

xenon detector had 3° circular collimation and an energy range of 1.5 - 35 keV 
while the argon detector had 2° x 8° collimation and an energy range of 
1.5 - 20 keV. 

Events were accepted as bona fide X-ray interactions if no coincidence 
between layers were present and if the collected charge were less than the 
upper threshold set to veto heavily ionizing events. The basic telemetry 
format consisted of contiguous 20.48 ms frames, each containing 64 words. of 
320 ps duration. The first event accepted during each word time was pulse-height, 
analyzed and digitized into one of 128 PHA channels. This event and 
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subsequent ones accepted during a word time -were also accumulated by a 
5-bit scaler with dead time of < 3 jusec per event. Rejected event counts 
were scaled every 20.48 ms and read out every frame. 

The attitude of the rocket was controlled by a STRAP IV (Shrewsberry , 
et al., 1973) pointing system, which, during the exposure to Cyg X-l, aimed 
the detectors to within 15 arcminutes of the programmed direction. This 
eliminated any appreciable modulation of the count rate due to motion of 
the rocket while pointing. 

III. DATA AND ANALYSIS 

Figure 1 displays a comparison of data from each of the three sources 
observed. Even though they are all considered collapsed objects in binary 
systems, they posses quite different X-ray temporal profiles. Her X-l 
exhibits clearly the 1.24 second pulsar period (Holt et al . , 1974), while 
Cyg X-3 displays no such obvious temporal structure. The X-ray profile 
of Cyg X-l shows considerable variability, qualitatively different from 
either of the other two sources. 

During the exposure of Cyg X-l the intensity experienced several 
periods of enhanced activity, especially from 318 to 319 seconds after 
launch (see Figure 2) . During this time the counting rate increased from 
the overall mean value of 1274 counts/second to 2188 counts/second averaged 
over 409.6 ms (20 data frames) encompassing the enhancement of about 1/3 
second duration. On three occasions during twenty milliseconds at the 
peak of the enhancement the instantaneous count rate for 1.28 ms samples 
increased to & 12 counts/1.28 ms (s 9375 counts/ second) (see Figure 3). 
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Poisson statistics, based on the mean rate of 2188 counts/second, would 

predict .01 such events during the 409,6 ms interval considered here; 

Using the enhancement and bursts as a guide, the entire Cyg X-l 

exposure was divided into 120 contiguous intervals of 409.6 ms duration, 

and each interval was divided into 320 1.28 ms bins. For each interval 

the mean count per bin was calculated and the counts for each bin noted. 

The expectation value N(n) for the number of bins containing n counts was 

determined for each interval using the mean count per bin for each interval 

and Poisson statistics. In each interval a count n* was picked such -that the 

S- 

expectation value for the number of bursts N* ^ -01, where n* was 

the lowest integer possible. All bins in the interval were then examined 

to see how many had n £ n*. The result of this analysis, summed over all 

120 intervals, was that in six intervals there existed at least one burst 

for which n > n*. These six intervals Contained eight bursts whereas 

the expectation value for the entire exposure is 0.73. The eight. bursts in 

. ' -4 

49.2 seconds of exposure indicate a duty cycle of 2 x 10 for detectable 
bursts. The eight bursts described above yielded a total of 85 counts of 
which 22 were pulse height analyzed over the band 2-10 keV (where detec- 
tion efficiency is high) . The efficiency-corrected average photon energy 
in this interval was 4.1 keV, indicating that the bursts may have a harder 
spectrum (but not one substantially softer) than the overall source emission. 

Adopting a technique similar to one brought to our attention by Oda 
(Oda, M., Matsuoka, M., Miyamoto, S., Ogawara, Y . , and Takagishi, K., 
Memorandum 1973) an attempt was made to estimate the burst duration of 
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the eight bursts described above assuming they had a simple "on-off" 
pulse shape with fixed width. Data from 6.72 ms encompassing each burst 
were used to evaluate the non-burst background in each case. The ratio 
of the net counts in the peak bin to the net counts in the adjacent 
spillover bin, superposed over all eight bursts, is consistent with a 
burst width of .9 ms for trial bins of both .96 ms and 1.28 ms. The 
same analysis of just the three bursts shown in Figure 3 was consistent 
with the same result. 

In evaluating the relative sensitivity of different rocket-borne 
experiments to millisecond bursts from Cyg X-l, it is important to note 
that if our detector area had been half as much as used in' this experiment 
the analysis employed here for detecting bursts would have recognized 
none of the bursts actually observed. 

In order to confirm the reliability of the above analysis, the data 
for Cyg X-3, (which had a mean count rate and exposure comparable to Cyg X-l) 
were examined for bursts and none were found. The rate of rejected (i.e., 
non X-ray) events was then examined during the time of the enhancement, 
and no significant variations were seen. This then minimizes the chance 
that the bursts were due to pick-up in the electronics from such things as 
attitude control valves firing or camera relays changing state. Finally, the 
combined rate from the inner layers of the argon detector was compared with 
the first layer rates from both detectors, and all three sets of data were 
similar within the statistical accuracy available. This eliminates the 
possibility that precipitating energetic electrons entering through the 
windows were responsible for the observed bursts. 



IV. DISCUSSION 


When there is disk accretion surrounding a black hole, fluctuations 
of the observed luminosity arising from turbulence may be expected 
(Shakura and Sunyaev, 1973) at a level given by 


AL „ M 2 

r 5 cr “ 


( 1 ) 


where M is the flux of matter, M is the flux corresponding to the 

cr 

Eddington critical luminosity, and a is an efficiency parameter (a ^ 1) 


for the mechanism of angular momentum transport. For a ~ 1 and M ~ M^, 
we note that AL is on the order of unity. The typical time of these 


L 

fluctuations is given by 

At ~ ) 3 ^ 2 seconds (2) 

f ot % 6GM 

where M is the mass of the black hole and R is the radial distance of the 
disturbance from the center of the black hole. For I£>R 0 = — ^ (the 
smallest stable orbit for the Schwarzchild metric), lOtfe, and 

o^l, we note that Atf > 10~ 3 seconds . 

The loss of energy (U) corresponding to a mass (M) reaching the 
innermost stable orbit is given by 

U = TIMc 2 (3) 

where!] is the efficiency of gravitational energy release (e.g.,^ .06 
for the case of the Schwarzschild metric and 11^.42 for the case of an 
extreme Kerr black hole) • The corresponding release of energy in X-radia- 
tion would be given by 

AL Atj = eU = eT]Mc 2 (4) 

where e is the efficiency of the X-radiatiori process relative to all energy- 
loss mechanisms. For thermal emission e approaches unity, whereas for non- 
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thermal emission e a 10 

The mass associated with the disturbances observed may be estimated 
by equation (4), viz: 

10 5 > °- 06Mc . 2 -^ 1 (5) 

ALfAtf 

where ALf is the luminosity increase during a fluctuation of duration Atf, 
and where we have taken n = .06. The lower limit corresponds to the situation 
where all energy transfer mechanisms channel their release into thermal X 
radiation (I.e., e ~ 1) . The upper limit corresponds to the non-thermal 
bremsstrahlung produced by particles transferring most of their energy to the 
ambient matter (i.e., e crl0~^) . For disturbances where At^ > 10”^ seconds 
(i.e. greater than the collision loss time), thermal radiation is likely to 
dominate, whereas for possible submillisecond structure, non-thermal radiation 
could be favored as the collision-loss time is approached. 

Using equation (5) we obtain that for the largest millisecond burst 
(containing 13 counts in a 1.28 ms bin compared with an expectation value 
of 2.8 counts) 

a 

7 x 10 19 o > M ^ 7 x 10 14 (6) 

whereas for the total associated time (~ 1/3 second) of enhanced activity, 
the corresponding composite mass involved becomes 

2 x 10 2i g > M Z 2 x lO^g. (7) 

These calculations assume the distance to Cyg X-l to be 2,5 kpc (Margon, 

Bowyer and Stone, 1973)* If we assume the infalling radial velocity is about 
10~ 2 times the Keplerian orbital velocity (see Pringle and Rees, 1972), then 
the 1/3 second duration of the observed enhancement corresponds to the infall 
time for matter from 4 R Q to R Q (the radius of the innermost stable orbit) 
about a 10 Mq black hole. The orbital velocities associated with this range 
of radii vary from B = .2 to B = .4. Particles undergoing velocity changes 
in this range are capable of producing the X-rays observed. 
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V. CONCLUSIONS 

We have evidence of the variability of Cyg X-l on time scales down 
to a millisecond, consistent with turbulence in disk accretion. The period 
of enhanced emission containing three of these bursts lasts for about 1/3 
second, which is consistent with the infall time for matter from 4 R Q to 
R 0 about a 10 M3 black hole. The changes in luminosity associated with 
these three bursts are greater than six times the mean luminosity of the 
entire exposure. When the entire exposure to Cyg X-l was examined, it 
was found that detectable bursts occurred ~ .02% of the time. 

Recent mass determinations of Cyg X-l (Brucato and Kristian, 1.973, 

Hutchings, et al., 1973, Bregman, et al., 1973) place its minimum mass 

greater than the 3.2 upper limit for neutron stars of any type (Leach 

and Ruffini, 1973). If one assumes that variations in emission come from 

turbulence in accretion, then millisecond bursts imply (see equation (2)) 

2 

that Cyg X-l must have a radius less than 10 km (e.g., less than the radius 
of a white dwarf) . This evidence combined with the mass determinations 
points toward a black hole nature for Cyg X-l. 
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FIGURE CAPTIONS 

Fifteen seconds of exposure to Her X-l , Cyg X-3, and Cyg X-l, 
respectively, binned every 20.48 ms. The inset shows the 
flight profile in counts 1 2 seconds versus seconds after launch. 
The horizontal arrows show where in the exposure to each source, 
the fifteen seconds of data occurred. 

The entire exposure to Cyg X-l as a function of time. The 
count rates are binned every 20.48 ms. Cyg X-l was acquired 
at 276 seconds after launch and the rocket doors started to 
close at 325 seconds after launch. The enhancement between 318 
and 319 seconds after launch contains three of the observed 
bursts . 

Eighty milliseconds of exposure to Cyg X-l containing the peak 
of the enhancement near 318 seconds after launch. The count 
rates are binned every 640 ji s. Bursts with ^ 12 counts/1.28 ms 


are shaded. 
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